Generation of large-scale vorticity in a homogeneous turbulence 

with a mean velocity shear 



Tov Elperin0 Nathan Kleeorin0 and Igor Rogachevski^ 

The Pearlstone Center for Aeronautical Engineering Studies, Department of Mechanical Engineering, 
Ben-Gunon University of the Negev, Beer-Sheva 84105, P. O. Box 653, Israel 

(Dated: February 2, 2008) 

An effect of a mean velocity sfiear on a turbulence and on the effective force which is determined 
by the gradient of Reynolds stresses is studied. Generation of a mean vorticity in a homogeneous 
incompressible turbulent flow with an imposed mean velocity shear due to an excitation of a large- 
scale instability is found. The instability is caused by a combined effect of the large-scale shear 
motions ("skew-induced" deflection of equilibrium mean vorticity) and "Reynolds stress- induced" 
generation of perturbations of mean vorticity. Spatial characteristics, such as the minimum size 
of the growing perturbations and the size of perturbations with the maximum growth rate, are 
determined. This instability and the dynamics of the mean vorticity are associated with the Prandtl's 
turbulent secondary flows. This instability is similar to the mean-field magnetic dynamo instability. 
Astrophysical applications of the obtained results are discussed. 

PACS numbers: 47.27.-1; 47.27.Nz 



I. INTRODUCTION 

Vorticity generation in turbulent and laminar flows was 
studied experimentally, theoretically and numerically in 
a number of publications (see, e.g., 0, 0, 0, S IE S 
i, H m in d ) • For instance, a mechanism of the 
vorticity production in laminar compressible fluid flows 
consists in the misalignment of pressure and density gra- 
dients 0, 01 . It was shown in that the vortic- 
ity generation represents a generic property of any slow 
nonadiabatic laminar gas flow. In incompressible flows 
this effect does not occur. The role of small-scale vor- 
ticity production in incompressible turbulent flows was 
discussed in [lo| . 

On the other hand, generation and dynamics of the 
mean vorticity are associated with turbulent secondary 
flows (see, e.g., P,ll,0,|^,fic«Z;i^J)- These flows, e.g., arise 
at the lateral boundaries of three-dimensional thin shear 
layers whereby longitudinal (streamwise) mean vorticity 
is generated by a lateral deflection or " skewing" of an ex- 
isting shear layer 0] . The skew-induced streamwise mean 
vorticity generation corresponds to Prandtl's first kind 
of secondary flows. In turbulent flows, e.g., in straight 
noncircular ducts, streamwise mean vorticity can be gen- 
erated by the Reynolds stresses. The latter is Prandtl's 
second kind of turbulent secondary flows, and it "has no 
counterpart in laminar flow and cannot be described by 
any turbulence model with an isotropic eddy viscocity" 

i 

In the present study we demonstrated that in a homo- 
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geneous incompressible turbulent flow with an imposed 
mean velocity shear a large-scale instability can be ex- 
cited which results in a mean vorticity production. This 
instability is caused by a combined effect of the large- 
scale shear motions (" skew- induced" deflection of equilib- 
rium mean vorticity) and " Reynolds stress- induced" gen- 
eration of perturbations of mean vorticity. The "skew- 

. . - is) . 

induced" deflection of equilibrium mean vorticity W is 
determined by (W'^'*^ • V)U-term in the equation for the 
mean vorticity, where U are perturbations of the mean 
velocity (see below) . The " Reynolds stress-induced" gen- 
eration of a mean vorticity is determined by V x where 
is an effective force caused by a gradient of Reynolds 
stresses. 

This instability is similar to the mean-field magnetic 
dynamo instability (see, e.g., 0|) which is caused by 
a combined effect of a nonuniform mean flow (differen- 
tial rotation or large-scale shear motions) and turbulence 
effects (helical turbulent motions which produce the a ef- 
fect or anisotropic turbulent motions which cause the 
"shear-current" effect [TH). 

This paper is organized as follows. In Section II the 
governing equations are formulated. In Section III the 
general form of the Reynolds stresses in a homogeneous 
turbulence with an imposed mean velocity shear is found 
using simple symmetry reasoning, and the mechanism for 
the large-scale instability caused by a combined effect 
of the large-scale shear motions and "Reynolds stress- 
induced" generation of perturbations of the mean vortic- 
ity is discussed. In Section IV the equation for the second 
moment of velocity fluctuations in a homogeneous tur- 
bulence with an imposed mean velocity shear is derived. 
This allows us to study an effect of a mean velocity shear 
on a turbulence and to calculate the effective force de- 
termined by the gradient of Reynolds stresses. Using the 
derived mean-field equation for vorticity we studied in 
Section IV the large-scale instability which causes the 
mean vorticity production. 



2 



II. THE GOVERNING EQUATIONS 

Our goal is to study an effect of mean velocity shear 
on a turbulence and on a dynamics of a mean vorticity. 
The equation for the evolution of vorticity W = Vxv 
reads 



aw 



Vx(vxW-i^VxW) , 



(1) 



where v is the fluid velocity with V • v = and f is 
the kinematic viscosity. This equation follows from the 
Navier-Stokes equation. In this study we use a mean 
field approach whereby the velocity and vorticity are sep- 
arated into the mean and fluctuating parts: v = U + u 
and W = W + w, the fluctuating parts have zero mean 
values, and U = (v) , W = (W) . Averaging Eq. over 
an ensemble of fluctuations we obtain an equation for the 
mean vorticity W 



dW 



Vx(UxW+ (uxw) -i/VxW) 



(2) 



Note that the effect of turbulence on the mean vorticity 
is determined by the Reynolds stresses (uiUj), because 



(3) 



and curl of the last term in Eq. Q vanishes. 

We consider a turbulent flow with an imposed mean ve- 
locity shear ViU*^'*\ where U^^-' is a steady state solution 
of the Navier-Stokes equation for the mean velocity field. 
In order to study a stability of this equilibrium we con- 
sider perturbations U of the mean velocity, i.e., the total 
mean velocity is U = U^'*-' -|-U. Similarly, the total mean 
vorticity is^W = W(") + W, where W^'*) = VxU(«) and 
W = VxU. Thus, the linearized equation for the small 
perturbations of the mean vorticity, W = W — W(*\ is 
given by 



aw 



= Vx(U(")xW-hUxW(") -h:F- ^VxW) , (4) 



where = — Vj[/ij(U) — fij (TJ^"')] is the effective force 
and fij — (uiUj). Equation Q is derived by subtracting 
Eq. Q written for W'^^ from the corresponding equation 
for the mean vorticity W. In order to obtain a closed 
system of equations in Section IV we derived an equation 
for the effective force Equation Q determines the 
dynamics of perturbations of the mean vorticity. In the 
next Sections we will show that under certain conditions 
the large-scale instability can be excited which causes the 
mean vorticity production. 



III. THE QUALITATIVE DESCRIPTION 

In this Section we discuss the mechanism of the 
large-scale instability. The mean velocity shear can 



LINEAR SHEAR FLOW = Sx KOLMOGOROV TURBULENCE u 



ANISOTROPIC "TANGLING" TURBULENCE oc - u ■ V U 



FIG. 1: Mechanism for tlie tangling turbulence: generation 
of anisotropic velocity fluctuations by tangling of the mean 
velocity gradients with the Kolmogorov-type turbulence. 



affect a turbulence. The reason is that additional 
strongly anisotropic velocity fluctuations can be gener- 
ated by tangling of the mean-velocity gradients with the 
Kolmogorov-type turbulence (see FIG. 1). The source 
of energy of this "tangling turbulence" is the energy of 
the Kolmogorov turbulence |l6|. The tangling turbu- 
lence is an universal phenomenon, e.g., it was introduced 
by Wheelon and Batchelor et al. for a pas- 
sive scalar and by Golitsyn and Moffatt |23| for a 
passive vector (magnetic field). Anisotropic fluctuations 
of a passive scalar (e.g., the number density of particles 
or temperature) are produced by tangling of gradients 
of the mean passive scalar field with a random veloc- 
ity field. Similarly, anisotropic magnetic fluctuations are 
generated by tangling of the mean magnetic field with 
the velocity fluctuations. The Reynolds stresses in a tur- 
bulent flow with a mean velocity shear is another exam- 
ple of a tangling turbulence. Indeed, they are strongly 
anisotropic in the presence of shear and have a steeper 
spectrur n (oc k~J^^ ) than a Kolmogorov turbulence (see, 
e.g., Tg", 21, 22, 23, "241). The anisotropic velocity fluctu- 
ations of tangling turbulence were studied first by Lumley 

The general form of the Reynolds stresses in a tur- 
bulent flow with a mean velocity shear can be obtained 
from simple symmetry reasoning. Indeed, the Reynolds 
stresses (uiUj) is a symmetric true tensor. In a tur- 
bulent flow with an imposed mean velocity shear, the 
Reynolds stresses depend on the true tensor VjUi, which 
can be written as a sum of the symmetric and anti- 
symmetric parts, i.e., WjUi = {dU).ij - {l/2)e^jkWk, 
where {dU)ij — {^iUj + V jJJi)/2 is the true tensor and 
W = VxU is the mean vorticity (pseudo- vector). We 
take into account the effect which is linear in perturba- 
tions {dU),j and W, where {dU),j = {ViUj + VjC/,)/2. 
Thus, the general form of the Reynolds stresses can be 
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found using the following true tensors: {dU)ij, Mij, iV^ 
Hij and Gij, where 







(5) 


% 


[^n-im (^^^ ^)mj ^rijmiS^U^ ^)m'i] 7 


(6) 






(7) 


Gij 




(8) 




= (V,;C7j^' + Vji7f''^)/2 and £,^7, is the fully 


an- 



tisymmetric Levi-Civita tensor (pseudo-tensor). There- 
fore, the Reynolds stresses have the following general 
form: 



+B2N,j+B3H,,+BiG,j), (9) 

where Bk are the unknown coefficients, Iq is the maxi- 
mum scale of turbulent motions, = IqUq/P is the tur- 
bulent viscosity with the factor /3 « 3 — 6, and uq is the 
characteristic turbulent velocity in the maximum scale 
of turbulent motions Iq. The parameter Iq in Eq. Q was 
introduced using dimensional arguments. The first term 
in RHS of Eq. describes the standard isotropic tur- 
bulent viscosity, whereas other terms are determined by 
fluctuations caused by the imposed velocity shear V^f/j . 

Let us study the evolution of the mean vorticity using 
Eqs. 0) and where Ti = — Vj/y (U) is the effective 
force. We consider a homogeneous divergence-free turbu- 
lence with a mean velocity shear, e.g., U^*^ — {0,Sx,0) 
and W*^"^ = (0, 0, S). For simplicity we use perturbations 
of the mean vorticity in the form W — {Wx{z), Wy{z), 0). 
Then Eq. Q can be written as 

dW 

= SWy+l^^W':, (10) 

^ = -asiiw':+^^w:^ , (11) 

where W'^ = d^W^^/dz^, a = (Si +2(^2 + B3) - 4B4)/4. 
In Eq. (fTn)l we took into account that l^Wy <C Wy, i.e., 
the characteristic scale Lw of the mean vorticity varia- 
tions is much larger than the maximum scale of turbulent 
motions Iq. This assumption corresponds to the mean- 
field approach. For derivation of Eqs. pU|) and (|11|) we 
used the identities presented in Appendix A. 

We seek for a solution of Eqs. pOfl and (|11|) in the form 
(X exp(7t -|- iifz). Thus, when a > perturbations of the 
mean vorticity can grow in time and the growth rate of 
the instability is given by 

-1 ^ ^/^SIqK -v^K"^ . (12) 

The maximum growth rate of perturbations of the mean 
vorticity, 7max = ot {S Iq)'^ / Aiy^ , is attained at K = Km — 
'JolSIqI'Iv^. The sufficient condition 7 > for the ex- 
citation of the instability reads Lw /lo > 27r/(/3 y^toS*), 



W^'^ = izS 




FIG. 2: Mechanism for "skew-induced" generation of pertur- 
bations of the mean vorticity Wx by quasi-inviscid deflection 
of the equilibrium mean vorticity W'''^ In particular, the 
mean vorticity Wx is generated by an interaction of pertur- 
bations of the mean vorticity Wy and the equilibrium mean 
vorticity wW, i.e., W^e^ cx (W'^' • V)l7,e, cx WyeyXW'-'K 

where Lyy = 2tt/K and we consider a weak velocity shear 
TqS < 1. 

Now let us discuss the mechanism of this instability 
using a terminology from The first term, SWy — 
(Wi'^) . V) Ux, in Eq. describes a "skew-induced" 
generation of perturbations of the mean vorticity Wx by 
quasi-inviscid defiection of the equilibrium mean vortic- 
ity W^''). In particular, the mean vorticity Wx is gen- 
erated from Wy by equilibrium shear motions with the 
mean vorticity W^*^ i.e., W^Gx oc (W^^^ • 'V)Uxex oc 
WyBx X W^"^ (see FIG. 2). Here e^, and Gy are the unit 
vectors along x and y axis. On the other hand, the first 
term, — aS'^g W^', in Eq. ((TT|l determines a "Reynolds 
stress-induced" generation of perturbations of the mean 
vorticity Wy by turbulent Reynolds stresses (see FIG. 
3). In particular, this term is determined by (Vx.^)y, 
where ^ is a gradient of Reynolds stresses. This implies 
that the mean vorticity Wy is generated by an effective 
anisotropic viscous term cx —Iq S W'x which is due to the 
equilibrium shear motions. The growth rate of this in- 
stability is caused by a combined effect of the sheared 
motions (" skew- induced" generation) and the "Reynolds 
stress-induced" generation of perturbations of the mean 
vorticity. 

This large-scale instability is similar to a mean-field 
magnetic dynamo instability (see, e.g., ^3)- Indeed, the 
first term in Eq. H10() is similar to the differential rota- 
tion (or large-scale shear motions) which causes a gen- 
eration of a toroidal mean magnetic field by a strech- 
ing of the poloidal mean magnetic field with the differ- 
ential rotation. On the other hand, the first term in 
Eq. Hll|l is similar to the a effect [lj|i or to the "shear- 
current" effect These effects result in generation of 
a poloidal mean magnetic field from the toroidal mean 
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FIG. 3: Mechanism for a "Reynolds stress-induced" genera- 
tion of perturbations of the mean vorticity Wy by turbulent 
Reynolds stresses. In particular, the mean vorticity Wy is gen- 
erated by an effective anisotropic viscous term oc —Iq S 
which is caused by the equilibrium shear motions. 



magnetic field. The a-effect is related with the hydrody- 
namic helicity of an inhomogeneous turbulent flow, while 
the "shear-current" effect occurs due to an interaction 
of the mean vorticity and electric current in a homoge- 
neous anisotropic turbulent flow of a conducting fluid. 
The magnetic dynamo instability is a combined effect of 
nonuniform mean flow (differential rotation or large-scale 
shear motions) and turbulence effects (helical turbulent 
motions which produce the a effect or anisotropic turbu- 
lent motions which cause "shear-current" effect). 

On the other hand, the magnetic dynamo instability 
is different from the instability of the mean vorticity 
although they are governed by similar equations. The 
mean vorticity W = V x U is directly determined by the 
velocity field U, while the magnetic field depends on the 
velocity field through the induction and Navier-Stokes 
equations. 



IV. EFFECT OF A MEAN VELOCITY SHEAR 
ON A TURBULENCE AND LARGE-SCALE 
INSTABILITY 

In this section we study quantitatively an effect of a 
mean velocity shear on a turbulence. This allows us to 
derive an equation for the effective force ^ and to study 
the dynamics of the mean vorticity. 



A. Method of derivations 

To study an effect of a mean velocity shear on a tur- 
bulence we used equation for fluctuations u(i, r) which is 



obtained by subtracting equation for the mean field from 
the corresponding equation for the total field: 



^.-,O.V,u-,u.V,U-^ 



+F(''*) -I- u 



N 



(13) 



where p are the pressure fluctuations, p is the fluid den- 
sity, F'*^*) is an external stirring force with a zero mean 
value, and Un — ((u • V)u) — (u • V)u -I- i^Au. We 
consider a turbulent flow with large Reynolds numbers 
(Re = Iquq/v ^ 1). We assumed that there is a sep- 
aration of scales, i.e., the maximum scale of turbulent 
motions Iq is much smaller then the characteristic scale 
of inhomogeneities of the mean fields. Using Eq. H13|) 
we derived equation for the second moment of turbu- 
lent velocity field /,j(k,R) = /(uj(k + K/2)uj(-k + 
K/2))exp(iK-R) dK: 



^/^J-(k,R) 
dt 



(N) 



(14) 



(see Appendix B), where 



(U) 



d 



^iq^jn^mp ~f~ ^im^jn^q | ^ pUq , (15) 



and R and K correspond to the large scales, and 
r and k to the small scales (see Appendix B), 

fcy = kikj/k'^, fj^^\'k,'R,) is the third moment ap- 
pearing due to the nonlinear term, V = 9/9R, 
Fy(k,R) = (#,(k,R)u,(-k,R)) + (u,(k,R)#,(-k,R)) 
and F(k,R,i) = -kx (kxF("*)(k, R))/fc2. Equation 
is written in a frame moving with a local velocity U of 
the mean fiow. In Eqs. (|14|l and (|15|l we neglected small 
terms which are of the order of 0(|V^U|). Note that 
Eqs. (|14l) and H15|) do not contain terms proportional to 
OdV^UI). 

The total mean velocity is U = U'"' + U, where we 
considered a turbulent fiow with an imposed mean ve- 
locity shear VitJ'^*^ . Now let us introduce a background 
turbulence with zero gradients of the mean fiuid velocity 
V,;C/j — 0. The background turbulence is determined by 



equation dff^'' /dt = Fij + flj'^^ where the superscript 
(0) corresponds to the background turbulence, and we 
assumed that the tensor Fij(k, R), which is determined 
by a stirring force, is independent of the mean veloc- 



(7V0) 



ity. Equation for the deviations fij 
ground turbulence is given by 



from the back- 



d{f 



= [/(U(^)) + /(tj)]/ + /(^) - / 



dt 



(NO) 



(16) 



where we used the following notations: / = 
/W ^ /f)(k,R), /(^o) ^ /f°^(k,R) 
/(">(k,R), and /(U)/ = /,,,„„(U)/™„(k,R). 



/»,(k,R) 
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Equation H16() for the deviations of the second moments 
in k-space contains the deviations of the third moments 
and a problem of closing the equations for the higher mo- 
ments arises. Various approximate methods have been 
proposed for the solution of problems of this type (see, 
e.g., [i^li^ll^l)- The simplest procedure is the r approx- 
imation which was widely used for study of different prob- 
lems of turbulent transport (see, e.g., [25, 28, 29, _30J). 
One of the simplest procedures, that allows us to express 
the deviations of the third moments f'^^^ — /(^"^^ in k- 
space in terms of that for the second moments f — f^^\ 
reads 

fiN)_fiNO) ^ -l-i^, (17) 
T{k) 

where T{k) is the correlation time of the turbulent veloc- 
ity field. Here we assumed that the time T(fc) is indepen- 
dent of the gradients of the mean fluid velocity because 
in the framework of the mean-field approach we may only 
consider a weak shear: ro|VC/| <^ 1, where tq — Iq/uq. 

The T-approximation is in general similar to Eddy 
Damped Quasi Normal Markowian (EDQNM) approx- 
imation. However there is a principle difference between 
these two approaches (see 25, 27] ). The EDQNM clo- 
sures do not relax to the equilibrium, and do not describe 
properly the motions in the equilibrium state. Within the 
EDQNM theory, there is no dynamically determined re- 
laxation time, and no slightly perturbed steady state can 
be approached |2^. In the r-approximation, the relax- 
ation time for small departures from equilibrium is deter- 
mined by the random motions in the equilibrium state, 
but not by the departure from equilibrium 25] . Analysis 
performed in showed that the r-approximation de- 
scribes the relaxation to the equilibrium state (the back- 
ground turbulence) more accurately than the EDQNM 
approach. 

Note that we applied the r-approximation l|17|) only 
to study the deviations from the background turbulence 
which are caused by the spatial derivatives of the mean 
velocity. The background turbulence is assumed to be 
known. Here we used the following model of the back- 
ground isotropic and homogeneous turbulence: 

f^PiKR) = -^/^,^.(k)f (fc) , (18) 

where -Pij(k) — Sij — kij, Sij is the Kronecker tensor, 
r(fc) ^ 2rof(fc), £{k) = -df{k)/dk, f{k) = (fc/fco)^"", 
1 < g < 3 is the exponent of the kinetic energy spectrum 
(e.g., q — 5/3 for Kolmogorov spectrum), and fco — l/^o- 



B. Equation for the second moment of velocity 
fluctuations 

We assume that the characteristic time of variation of 
the second moment /ij(k, R) is substantially larger than 



the correlation time r(fc) for all turbulence scales. Thus 
in a steady-state Eq. reads 

L{I - = r{k) + /(U)] /(O) , (19) 

where L = Lym„ = SimSjn - ^(fc) [Jy„in(U*^'''^) + 
-fiimn(U)], and we used Eq. l(T7|) . The solution of Eq. itT^ 
yields the second moment / = /y(k, R): 

/ « /(U(^)) + T{k) [/(U) + /(U(^)) r(fc) /(U) 

+/(U)r(fc)/(U(^))]/W , (20) 

where /(U^"') ^ /(") + r(fc)[/(U(")) + 
/(U(''))r(fc)/(U(''))]/(o). In Eq. ^ we neglected 
terms which are of the order of 0(]VU]^) and 
0( jVU^^'l^jVU]). The first term in the equation 
for /(U^"^) is independent of the mean velocity shear 
and it describes the background turbulence, while the 
second and the third terms in this equation determine 
an effect of the mean velocity shear on turbulence. 

C. Effective force 

Equation (|2()(l allows us to determine the effective 
force: T., = -V^ / (k, R) dk, where / = / - /(U^^)), 
and we used notation / = /,y (k, R). The integration in k- 
space yields the second moment f ij (R) = J fij (k, R) dk: 

^(R) = -2v^{dU\j-ll[ACiM,, 

+C2{N,,+H,j)+C^G,j], (21) 

where Ci = 8(g2 - I3q -|- 40)/315, Ca = 2(6 - 7g)/45, 
C3 = —2{q -I- 2)/45, and the tensors Afy , iV^ , Hij and 
Gij are determined by Eqs. (0)-®- In Eq. 1^1)1 we omit- 
ted terms cx Sij because they do not contribute to Vx:F 
(see Eq. Q for perturbations of the mean vorticity W) . 
To derive Eq. iglj we used the identities presented in 
Appendix A. Equations © and if^ yield Bi = 4Ci, 
52 = ^3 = C2 and B4 = C3. 

Note that the mean velocity gradient ViU*^*-' causes 
generation of anisotropic velocity fluctuations (tangling 
turbulence). Inhomogeneities of perturbations of the 
mean velocity U produce additional velocity fluctuations, 
so that the Reynolds stresses /ij(R) are the result of a 
combined effect of two types of velocity fluctuations pro- 
duced by the tangling of mean gradients ViU'^*-' and ViU 
by a small-scale Kolmogorov turbulence. Equation l|21l) 
allows to determine the effective force J^i = — Vj/ij(R). 

D. The large-scale instability in a homogeneous 
turbulence with a mean velocity shear 

Let us study the evolution of the mean vorticity us- 
ing Eq. H21|) for the Reynolds stresses. Consider a ho- 
mogeneous turbulence with a mean velocity shear, e.g.. 
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UW = (0,S'a;,0) and W(^) = (0,0,5). For simplic- 
ity we consider perturbations of the mean vorticity in 
the form W = {W^{z),Wy{z),0). Then Eq. (jH) reduces 
to Eqs. lUni and (jTTl), where a = Ci + C2 - C3 = 
4(2g2 _ 47? + 108)/315, and we used Eq. (EU- We 
seek for a solution of Eqs. H1Q|1 and in the form 
cx exp(7t + iKz). Thus, the growth rate of perturbations 
of the mean vorticity is given by 7 = ^/aSlQ K — . 
The maximum growth rate of perturbations of the mean 
vorticity, 7max — oi{SIqY /4v^ 0.1 (3 S'^tq, is attained 
at K = = y/aSlol2v^. Here we used that for 

a Kolmogorov spectrum {q = 5/3) of the background 
turbulence, the factor a « 0.45. The sufficient condi- 
tion 7 > for the excitation of the instability reads 
Lwlh > ^-K / {P ^/a tqS) . Since tqS <^ 1 (we considered 
a weak velocity shear) , the scale Lw ^ ^0 and, therefore, 
there is a separation of scales. The mechanism of this 
instability is discussed in Section III and is associated 
with a combined effect of the "skew-induced" deflection 
of equilibrium mean vorticity due to the sheared motions 
and the "Reynolds stress-induced" generation of pertur- 
bations of mean vorticity. 



V. CONCLUSIONS AND APPLICATIONS 

We discussed an effect of a mean velocity shear on 
a turbulence and on the effective force which is deter- 
mined by the gradient of Reynolds stresses. We demon- 
strated that in a homogeneous incompressible turbulent 
flow with an imposed mean velocity shear a large-scale 
instability can be excited which results in a mean vortic- 
ity production. This instability is caused by a combined 
effect of the large-scale shear motions ("skew-induced" 
deflection of equilibrium mean vorticity) and "Reynolds 
stress-induced" generation of perturbations of mean vor- 
ticity. We determined the spatial characteristics, such as 
the minimum size of the growing perturbations and the 
size of perturbations with the maximum growth rate. 

The analyzed effect of the mean vorticity produc- 
tion may be of relevance in different turbulent indus- 
trial, environmental and astrophysical flows (see, e.g., 
Slil5ll6ll7l|8l|31l|3aiMil). Thus, e.g., the sug- 
gested mechanism can be used in the analysis of the 
flows associated with the Prandtl's turbulent secondary 
flows (see, e.g., |lllllll[lll). These flows, e.g., arise 
in straight noncircular ducts, at the lateral boundaries 
of three-dimensional thin shear layers, etc. The simple 
model considered in the present paper can mimic the 
flows associated with turbulent secondary flows. 

The obtained results may be also important in astro- 
physics, e.g., in extragalactic clusters and in interstellar 
clouds. The extragalactic clusters are nonrotating ob- 
jects with a homogeneous turbulence in the center of a 
extragalactic cluster. Sheared motions between interact- 
ing clusters can cause an excitation of the large-scale in- 
stability which results in the mean vorticity production 
and formation of large-scale vortices. Dust particles can 



be trapped by these vortices to enhance agglomeration 
of material and formation of particles inhomogeneities 
H^HEI^' The sheared motions can also occur between 
interacting interstellar clouds, whereby the dynamics of 
the mean vorticity is important. 
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APPENDIX A: IDENTITIES USED FOR 
DERIVATION OF EQS. ^0^, AND (I2T]) 

To derive Eqs. H10|) and we used the following 
identities: 

VxM = {SKyA){Wy,W^,0) , 

Vx(UxW(")) = S{Wy,-W^,0), 

Vx(U('*)xW) = S{0,W^,0) , 

VxJ = {SKy2){Wy,W^,0) , 

Vx[(W(^) • V)W] - SK^{Wy,-W^,0) , 

where M, = V^My, J^ = Vj [W^„(£„„„(ai7(^))„y + 
£njmidU'''^^)m.i)], and we also took into account that 

V,[W^''\emm{dU)mj + (<9C/)™)] 
= [Vi(W(') • W) - (W^'') • V)VK,]/2 . 

To derive Eq. (|21|l we used the following identities for 
the integration over the angles in k-space: 



T- 
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-J 


^ijmnpq dCl 








-\- ^ijtnn^pq 



47r 
105' 



^rnnpq^ij 



— 47r 

Tym„pg(V™C7^"))(VpC/g) = — (4Mij +%Mpp) , 



kijmn^pq ^-Hd. ^ijfjin — ^ij^mn ~^ ^im^nj H" ^in^mj • 



APPENDIX B: DERIVATION OF EQ. (ITU) 

In order to derive Eq. (|I4|I we use a two-scale approach, 
i.e., a correlation function is written as follows 

(ui(x)Mj(y)) = y"(uj(ki)uj(k2))exp[i(krx 



7 



+k2-y)] dki dk2 
= y /ij (k,R)exp(ik-r)fik , 

/,,(k,R) = J{u,{k + K/2)uj{-k + K/2)) 
X exp(iK-R) dK 

(see, e.g., [s^ H^), where R and K correspond to the 
large scales, and r and k to the small scales, i.e., R = 
(x + y)/2, r = x-y, K = kj + kj, k = (ki - 
k2)/2. This implies that we assumed that there exists a 
separation of scales, i.e., the maximum scale of turbulent 
motions lo is much smaller then the characteristic scale 
of inhomogeneities of the mean fields. 
Now we calculate 

. (P„(k,)^^U,(k2)) 

+ iuMPM^^^) , (Bl) 



dt 



where we multiplied equation of motion rewritten 
in k-space by Pij{k) = 5ij — kij in order to exclude the 
pressure term from the equation of motion, Sij is the 
Kronecker tensor and kij = kikj/k'^. Thus, the equation 
for fij(k, R) is given by Eq. ((Tl|l . 

For the derivation of Eq. (|14|) we used the following 
identity 



ik^ I /y (k~ -Q,K- Q)C/p(Q)exp(iK-R)dKdQ 



1 
2 



|(V,V,i7p) . 



dks 



4 V dk. 



(B2) 



To derive Eq. ljB2|) we multiply the equation 
V • u = 0, written in k-space for Mi(ki — Q), by 
Mj(k2)[/p(Q) exp(iK-R), and integrate over K and Q, 
and average over ensemble of velocity fluctuations. Here 
ki = k + K/2 and ka = -k + K/2. This yields 

^i(k, + hi, - (u,(k + - Q)ii,(-k + Ik)) 
xC7p(Q)exp(iK-R)dKdQ = . (B3) 

Next, we introduce new variables: ki = k + K/2 — Q, 
k2 = -k + K/2 and k = (ki - k2)/2 = k - Q/2, K = 
ki + k2 = K — Q. This allows us to rewrite Eq. I|B3|) in 
the form 



» ( fc, + -K, - ] /„ (k - ^Q, K - Q)C/p(Q) 



X exp(iK-R)dKdQ = . (B4) 
Since |Q| <C |k| we can use the Taylor expansion 



/,,(k - Q/2, K - Q) ~ /y(k, K - Q) 



ia/y(k,K-Q) 



2 dks 
We also use the following identities: 



Qs + 0{Q^ 



(B5) 



[/,,(k,R)f/p(R)]K = J /,,(k,K - Q)C/p(Q)dQ , 

Vp[/y(k,R);7p(R)] = J zifp[/y (k,R)f7p(R)]K 
X exp (iK-R) dK . (B6) 

Therefore, Eqs. (|B4l)-||B6ll yield Eq. l|52|l . 
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